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Abstract 

We consider a bivariate Cramer-Lundberg-type risk reserve process with the special feature that each in¬ 
surance company agrees to cover the deficit of the other. It is assumed that the capital transfers between 
the companies are instantaneous and incur a certain proportional cost, and that ruin occurs when neither 
company can cover the deficit of the other. We study the survival probability as a function of initial capitals 
and express its bivariate transform through two univariate boundary transforms, where one of the initial 
capitals is fixed at 0. We identify these boundary transforms in the case when claims arriving at each 
company form two independent processes. The expressions are in terms of Wiener-Hopf factors associated 
to two auxiliary compound Poisson processes. The case of non-mutual agreement is also considered. The 
proposed model shares some features of a contingent surplus note instrument and may be of interest in the 
context of crisis management. 

Keywords: Two-dimensional risk model, survival probability, coupled processor model, Wiener-Hopf 
factorization, surplus note, mutual insurance 


1. Introduction 

Insurance companies cannot operate in isolation from financial markets or from other insurance and rein¬ 
surance companies. Hence it is important to understand the effect of interaction on the main characteristics 
of an insurance company. Multivariate risk models, however, present a serious mathematical challenge with 
few explicit results up to date, see [1, Ch. XHI.9]. This paper focuses on a nonstandard, but rather general 
bivariate risk model and provides an exact analytic study of the corresponding survival probability, borrow¬ 
ing some ideas from the analysis of a somewhat related queueing problem in [7] . The proposed model may 
be of interest in the context of crisis management due to its relation to contingent surplus notes and mutual 
insurance. Moreover, it allows for an explicit structural result without imposing overly strict assumptions 
such as proportional or dominating claims, see Section 1.1. 

We consider a bivariate Cramer-Lundberg risk process as a model of surplus of two insurance companies 
(or two lines of one insurance business). The special feature of our model is that the companies have a 
mutual agreement to cover the deficit of each other. More precisely, if company 1 gets ruined, with its 
capital decreasing to a value —x < 0, then company 2 compensates this deficit, bringing the capital of 
company 1 back to 0. However, this comes at a price; a unit of capital received by company 1 requires 
ri > 1 from company 2 (cf. Figure 1). If this would cause the capital of company 2 to go below 0, then both 
companies are said to be ruined. Similarly, if company 2 gets a deficit —y < 0, company 1 compensates this 
deficit, but its capital reduces by r 2 y > y, and if this would cause the capital of company 1 to go below 0, 
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then again both companies are said to be ruined. Finally, ruin may also be caused by a single event bringing 
the surplus processes of both companies below 0. It may be more realistic to assume that if one company can 
not save the other from ruin then it does not transfer any capital at that instant and continues to operate. 
Note, however, that survival of both companies in this set-up corresponds to our previous notion of survival. 



Figure 1: Illustration of mutual deficit coverage (no common shocks). 

Our main goal in this study is to provide an exact analysis of the (Laplace transform of the) probability 
of survival (i.e., ruin never occurs) as a function of the vector {u,v) of initial capitals. We do this 

by (i) expressing the two-dimensional Laplace transform F(si,S 2 ) of (f>{u,v) in terms of the transforms of 
(/)(u, 0) and (/)(0, n), and (ii) determining the latter two transforms by solving a Wiener-Hopf boundary value 
problem in the case of independent claim streams. In the latter step, a key role is played by the Wiener-Hopf 
factorization of two auxiliary compound Poisson processes. 

In our terminology, if one company cannot save the other then both are declared ruined, whereas in 
practice it is more likely that the company with a positive capital pays out all it has, but then continues its 
operation. Hence it may be more appropriate to call (j) the probability of survival of both companies, or the 
probability that external help is never needed, see Remark 3. One may also notice that it is clearly better 
(with respect to survival) to merge the two lines eliminating transaction costs. There may be cases, however, 
where a merger is not possible due to legal, regulatory or other issues. Outside of an insurance context, the 
two lines may be two separate physical entities such as water reservoirs or energy sources. Furthermore, 
one may consider the case where < 1 for at least one company, see Remark 2, so that merging may not 
be optimal. This may correspond to the case where part of the deficit is written off or is covered by some 
other fund. 

Our model resembles two insurance companies with a mutual insurance fund, which is used to cover the 
deficit of each of them. In our case, however, there is no separate fund - it is the other company which 
provides the capital to cover the deficit. Another related notion is that of a contingent surplus note which 
is a form of a CoCo bond issued by an insurance company. This bond pays a higher coupon, because of 
the risk that it is converted into surplus if a trigger event (e.g. deficit) occurs. Finally, we may consider 
an insurance company and a governmental fund which is used to save companies from default. In this case 
capital transfers are only possible from the fund to the company; this particular scenario is discussed in 
Section 7. 

The above-sketched risk model bears some resemblance to a two-dimensional queueing model of two 
coupled processors. This model features two M/G/l queues, each of which, in isolation, is known to be 
the dual of a Cramer-Lundberg insurance risk model, cf. [1, Ch. L4a]. Just like in our risk model, the two 
processors are coupled by the agreement to help each other. When one of the two M/G/l queues becomes 
empty, the service speed of the other server - say, server i - increases from 1 to > 1, i.e. server i is 
being helped by the idle server. It should be stressed that this similarity is rather loose and that there is 
no clear duality relation between these bivariate risk and queueing models. Strikingly, the crucial ideas of 
the analysis of the coupled processor model in [7] apply to our present setup as well. Moreover, solutions 
to both problems are based on the same Wiener-Hopf factors, which hints that there might be a certain 
duality relation between the two. 
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1.1. Related literature 

Despite their obvious relevance, exact analytic studies of multidimensional risk reserve processes are 
scarce in the insurance literature. A special, important case is the setting of proportional reinsurance, which 
was studied in Avram et al. [3]. There it is assumed that there is a single arrival process, and the claims 
are proportionally split among two reserves. In this case, the two-dimensional exit (ruin) problem becomes 
a one-dimensional first-passage problem above a piece-wise linear barrier. Badescu et al. [4] have extended 
this model by allowing a dedicated arrival stream of claims into only one of the insurance lines. They show 
that the transform of the time to ruin of at least one of the reserve processes can be derived using similar 
ideas as in [3]. 

An early attempt to assess multivariate risk measures can be found in Sundt [22], where multivariate 
Panjer recursions are developed which are then used to compute the distribution of the aggregate claim 
process, assuming simultaneous claim events and discrete claim sizes. Other approaches are deriving integro- 
differential equations for the various measures of risk and then iterating these equations to find numerical 
approximations [10, 16], or computing bounds for the different types of ruin probabilities that can occur in 
a setting where more than one insurance line is considered, see [8, 9]. In an attempt to solve the integro- 
differential equations that arise from such models, Chan et al. [10] derive a Riemann-Hilbert boundary value 
problem for the bivariate Laplace transform of the joint survival function (see [5] for details about such 
problems arising in the context of risk and queueing theory and the book [11] for an extended analysis of 
similar models in queueing). However, this functional equation for the Laplace transform is not solved in 
[10]. In [5] a similar functional equation is taken as a departure point, and it is explained how one can 
find transforms of ruin related performance measures via solutions of the above mentioned boundary value 
problems. It is also shown that the boundary value problem has an explicit solution in terms of transforms, 
if the claim sizes are ordered. In [6] this is generalized to the case in which the claim amounts are also 
correlated with the time elapsed since the previous claim arrival. 

Bivariate models where one company can transfer its capital to the other have already been considered 
in the literature. Recently, Avram and Pistorius [2] proposed a model of an insurance company which splits 
its premiums between a reinsurance/investement fund and a reserves fund necessary for paying claims. In 
their setting only the second fund receives claims, and hence all capital transfers are one way: from the first 
fund to the second. Another example is a capital-exchange agreement from [20, Ch. 4], where two insurers 
pay dividends according to a barrier strategy and the dividends of one insurer are transferred to the other 
unless the other is also fully capitalized. This work resulted in systems of integro-differential equations for 
the expected time of ruin and expected discounted dividends, which are hard to solve even in the case of 
exponential claims. 

Finally, we briefly list related contributions in the queueing context. The joint queue length distribution 
of the coupled processor model has been derived by Fayolle and lasnogorodski [15], in the case that the 
service time distributions at both queues are exponential. In their pioneering paper, they showed how the 
generating function of the joint steady-state queue length distribution can be obtained by solving a Riemann- 
Hilbert boundary value problem. Cohen and Boxma [11] generalized this queueing model by allowing general 
service time distributions. They obtained the Laplace-Stieltjes transform of the joint steady-state workload 
distribution by solving a Wiener-Hopf boundary value problem. In [7] the model of [11] was extended by 
considering a pair of coupled queues driven by independent spectrally positive Levy processes and a compact 
solution was obtained. There the model was also linked to a two-server fluid network. 

1.2. Organization of the paper 

In Section 2 we describe the model in detail. In Section 3 we derive an integral equation for the survival 
probability as a function of the vector {u,v) of initial amounts of capital. Section 4 is devoted 

to the derivation of a so-called kernel equation for the two-dimensional Laplace transform of (j){u,v), see 
Proposition 1. After a brief discussion of the net profit condition, in Section 5, we solve the kernel equation 
in the case of independent claim streams in Section 6. The main result is given in Theorem 1, which is then 
illustrated by two simple examples. In Section 7 we specialize to a model of non-mutual insurance, where 
help is provided by only one company (or government). Some open problems are discussed in Section 8. 
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2. The model 


Consider a bivariate Cramer-Lundberg risk process {u + Xi(t),v + X 2 {t)),t > 0 with initial capitals 
M, u > 0, premium rates Ci > 0, claim arrival rate A and a joint claim size distribution /j,(dxi, dx 2 ). In other 
words, 

N{t) 

(Xiit), X2{t)) = {ci,C2)t— (Cl,i,C 2 ,i), 

i-1 

where {Ci^i,C 2 ,i),i S N are iid distributed according to fi, and N{t) is an independent Poisson process of 
rate A. Without loss of generality we can assume that /r does not have mass at (0, 0). Next, we define the 
corresponding bivariate Laplace exponent by 

ip{si,S2) := ( 1 ) 

= sici + S2C2 - A dxa)^ . 

Remark 1. Note that this model incorporates the case of two independent risk processes with claim arrival 
rates Xi and claim size distributions /ii(dx) with the relation 

A = Ai + Aa, dy) = yi{dx)6o{dy) + y 2 {dy)So{dx), 

Ai A2 Ai + A2 

where <5o denotes the Dirac point mass at 0. This yields ipisi, S 2 ) = V’i(si) + ^' 2 ( 52 ) with 


'tj}i{s) := logEe®^'^^^ = sci - A* 




z = 1,2. 


Throughout this work we assume that the companies have a mutual agreement of deficit coverage: when 
a claim brings the surplus of a company z € {1,2} below 0 then the other company j = 3 — z transfers just 
enough capital to bring it to 0 so that company z can continue its operation. If company j has not enough 
capital then ruin occurs. We also assume that a capital transfer incurs some proportional cost: a unit of 
capital received by i requires Vi > 1 capital from j. It is noted that in the case of no transaction costs our 
risk problem reduces to a classical one-dimensional problem obtained by considering the total surplus, see 
Section 6.2. 

Remark 2. All of the results below hold for any > 0 such that rir 2 > 1. Even the latter condition can 
be removed, except that then the second statement in Proposition 2 may no longer he true. 

It is convenient to consider a more general model defined for all times < > 0, where ruin is avoided using 
capital from external sources. More precisely, if company j has insufficient funds to save z then it transfers 
to z all its capital and the rest of the required capital is taken from the shareholders of company z (or some 
other external source). Mathematically, we can describe the bivariate surplus process (S'i(t), «S' 2 (t)) in the 
following way: 

Si{t) = u + Xi{t) — r 2 L 2 {t) -|- Li{t) + Ei(t), 

^ 2 ( 1 ) = u -|- X2{t) — riLi{t) -|- L2{t) + i?2(t), (2) 

where Li(f) and Ei{t) are the cumulative amounts of capital received by the z-th company from the other 
company and from an external source respectively; Li(0) = Ei(0) = 0. Note that the time of ruin in the 
original model is given by 

T = infjt >0:Ei{t)+ E 2 {t) > 0}. (3) 

In this work we only address the probability of survival 4>{u, v) := Pu.iiCt = 00 ). Note that in our setup the 
concepts of ruin and survival have unambiguous meaning, which should be compared to various possible 
ruin concepts in standard bivariate risk models, see [10]. In the following we provide two mathematical 
definitions of the processes involved in (2): A recursive definition and a definition as a modified Skorokhod 
problem. 
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Remark 3. We note that the external sources Ei are introduced only for modeling convenience. Throughout 
this work survival means survival in the original model which corresponds to never using external capital 
in the extended model. Hence Ei never appear in the following apart from the discussion of the net profit 
condition in Section 5 and direct calculations in Section 6.2. Nevertheless, we believe that it is important to 
have an explicit model definition as given by (2) and (3), see also Section 8. 

2.1. Recursive definition 

Let us remark that the processes Li, Ei are piece-wise constant and non-decreasing, and can have jumps 
at claim arrival times ti,t 2 ,... only. One may define these processes recursively by considering their jumps 
at time for n = 1,2,... More precisely, with Xi = Si{tn—) + {Xi{tn) — Xi{tn—)) we have the following 
cases: 


• Xi,X 2 > 0: Li,Ei have no jumps (at t„), 

• a;i < 0,a:2 -I- rixi > 0: Li has a jump —xi (and no other jumps at t„), 

• a ;2 < 0, a:i -I- r 2 X 2 > 0: L 2 has a jump —X 2 , 

• xi < 0, X 2 > 0, X 2 + rixi < 0: Li has a jump X 2 /ri and Ei has a jump —xi — X 2 lr\, 

• 0:2 < 0, a:i > 0, xi -I- r 2 X 2 < 0: L 2 has a jump Xi/r 2 and E 2 has a jump —X 2 — xxjr 2 , 

• xi < 0 ,X 2 < 0: Ei has a jump —Xi for * = 1, 2. 

2.2. Modified Skorokhod problem definition 

The following equivalent specification of our model closely resembles a multidimensional Skorokhod 
problem, see [17, 18] building upon Skorokhod’s one dimensional construction [21]. It is noted that for the 
classical bivariate Skorokhod problem one must assume that rir 2 < 1 which is clearly violated in our setup. 
We consider the equations in (2) with the following additional requirements: 

• Si > 0 and Li, Ei are piece-wise constant and non-decreasing, 

• the jump times of Li are contained in {t > 0 : Si(t) = 0}, 

• the jump times of Ei are contained in {t > 0 : Si{t) = S 2 {t) = 0}, 

• out of Li,L 2 , Ei,E 2 the only pairs which can jump simultaneously are {Li,Ei), {L 2 , E 2 ) and {Ei, E 2 ). 

It is easy to check that these requirements are satisfied by the recursive construction from Section 2.1, and 
that no other choices of Li,Ei and thus Si satisfy these conditions. 

A few words about the above requirements. The second states that a company receives capital when in 
need, and the third states that external capital can be used only when both companies are in need (or one is 
at 0 after trying to save the other). The final condition is necessary to prevent redundant capital transfers. 

3. The survival probability 

As stated in Remark 3, we now ignore the possibility of funding by external sources. Consider (/)(u, v) = 
IPii,j;(t = 00 ), the probability of survival for the initial capitals u,v>0, and put for convenience 4>{u, u) = 0 
if M < 0 or u < 0. It is immediate from sample path comparison that <fi{u, v) is non-decreasing in both u, v. 
Furthermore, 4>{u, v) is a (Lipschitz) continuous function on which follows from u) > (1 — Xh)(f){u + 
Cih, V -I- C 2 / 1 ) -I- o{h) as h fO implying 

4>{u + Cl ft., V + C 2 h) — <f>{u, v) < \h + o{h), (4) 
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where o{h) does not depend on m , v. Moreover, our model definition implies the following equation 

v) = o{h) + (1 — Xh)(j){u + Cih, V + C 2 h) 

+ Xh jj (j){u — x,v — y)ii{dx,dy) 


( 5 ) 


x^[0,u\,y^[0,v\ 


+ 


+ 


Xh JJ (j){u - X + r 2 {v - y),0)fj,{dx,dy) 

[0,n] ,y>v 

Xh JJ (f>{0,v — y + ri{u — x))ij{dx,dy). 


x'>u,yG[0,v] 


In fact, this equation requires some careful considerations, and so its proof is given in the Appendix. 
Equation (5) implies that the directional derivative 


D(c^,c^)(t>{u,v) = hm 


(j){u + Ci/l, V + C2I1) — v) 


exists and the following equation holds true. 

(j){u,v) = + JJ (j){u- x,v -y)y{dx,dy) 


xG[0,'u],yG[0,'u] 


+ 


JJ (j){u-X + r 2 {v - y),0)n{dx,dy) 

xG [0,ii] ,y>'v 

JJ - y + ri{u-x))y,{dx,dy). 


x'>u,y^ 


Remark 4. If the partial derivatives (pi and p 2 exist and are continuous at (u,v) then D(^ci,c2)4'{xLtV) = 
ci(pi(u, v) + C 2 p 2 {u, v). This can not hold in general when fi has an atom at {u, v), as can be seen from (5). 


4. The kernel equation 


Equation (5) can be equivalently formulated in terms of Laplace transforms. This results in a very simple 
identity given by the following Proposition. 

Proposition 1. Let F{si, S 2 ) ■= fg 2 e ^^'"p(u,v)dudv and 


pOO 

Fiisi) = / e~''^'^(p{u,0)du, 

Jo 


pOO 

^ 2 ( 52 )= / e~^^'^(j){ 0 ,v)dv. 

Jo 


Then for all Si, S 2 >0 it holds that 




32 - r2Si 


Si - riS2 


( 6 ) 


This result identifies the transform of p{u, v) up to the unknown functions Fi{si). Equation (6) resembles 
the so-called basic adjoint relation in the semimartingale reflected Brownian motion (SRBM) literature, see 
e.g. [23, Eq. (3.2)], with the class of functions f{u,v) = queueing literature this type of 

functional equation is sometimes called ‘kernel equation’, see e.g. [11, Ch. III.3] and [7, Prop. Ij. 
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Proof of Proposition 1. We simply take the transform of ( 5 ) and let h 0. Let us first consider the transform 
of 4>(u + cih, V + C2h) — (p(u, v): 


e + cih,v + C2h)dudv — F{si, S2) 


= u)dud^; - F(si, S2) 

jmi 

^ (^^siCih+S2C2h _ S2) - [ l{u<cih or u)dMdu. 


By dominated convergence and continuity of (j) we have 


lim — 
/40 h , 


n >0 'J v<C2h 


e (l){u, v)dudv 


= C2 [ e ®^“lim^— f e ‘‘^'"(j){u,v)dudv = C2Fi{si). 
Ju >0 ^2/1 Jv<C2h 


Similarly, we have 


e ^^'"(p{u,v)dudv = CiF2(s2) 


^ — SiU—S2V 


lim - / / 

h^O fl JuKc\h J t»>0 

(j){u^v)dudv = o{h). Combining these we obtain 


fu<cih fv<C 2 h ' 

lim - / ((/)(u + cih, V + C2h) — (l){u, u))dudu = (siCi + S2C2)F{si, S2) — C2Fi{si) — ciJ^2(s2)- 

HO ti 

It is left to consider the transforms of the integrals in ( 5 ). Using Fubini’s theorem we get 

/ - a:,i; - j/)/r(da;,dj/)dwdt; 

J«.\ 

= / x,v - y)dudvfj,{dx,dy) 

JRI JRl 

= [ [ e~^^^^'^^^~^'^^'’~^^^(j){u,v)dudvfi{dx,dy) = G{si,S2)F{si,S2), 

JR^ JRl 

where G(si,S2) = /r 2 dy). 

Next, using the substitution u' = u — x + r2{v — y) we obtain 

/ l{x<u,y>v}(t>{u - X + r2{v - y), 0 )p{dx,dy)dudv 

JR'^. 


— S 1 U—S 2 V 


/ / ^{u>x,v<v}(i -a; + r2(?;-y),0)dudi;/r(da::,dy) 

Jr\ 

JRI JRl 

[ - (1 

Jr 2 S2 — Sir2 


P-(S2-Sl»'2)j/ 


)fj.{dx,dy) 


Fi{s,) 
S2 - Sir2 


(G(si,sir2 ) - G(si,S2 )), 


( 7 ) 


and a similar equation for the last integral in ( 5 ). 
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Finally, we combine all the terms to get 


0 — (siCi + S2C2)F{si, S 2 ) — 01 ^ 2 ( 52 ) ~ C 2 Fi{si) 

+ X (^—F{si, S 2 ) + G(si, S2)F{si, S 2 ) 

H- - — '^—{G{si, Sir2) — G{si, S 2 )) H-^S 2 ) — G{si, 52 ))^ ■ 

S2 - Sl?’2 Si - S2?'l / 

Some simple manipulations using (1) show that this equation coincides with (6) and so we are done. □ 

The following Corollary for two independent driving processes is immediate. 

Corollary 1 . If Xi and X2 are independent proeesses given in Remark 1 then 

( I ( \ ^ I ( \ V’2(s2)-^/’2(?’2Si) , . V'l(si) - V'l(?'lS2) p 

[Viisi) +i’2is2))F{si,S2) = - Fi(si) H-^ 2 ( 32 ) (8) 

S2 - r2Si Si - riS2 

for all si,S2 > 0 . 

5 . The net profit condition 

Let Hi = EXi(l) G [— 00 , 00 ) be the average drift of the driving process Xi. In the following we assume 
that the net profit condition holds 

Hi + T 2 H 2 >0, H2 + riHi > 0. (9) 

The importance of this condition is explained by the following result. 

Proposition 2 . The following dichotomy is true for ri G (0, 00 ) with rir2 > 1.' 

• if (9) holds then (/)(oo,0) = ^(0, 00 ) = 0(oo,oo) = 1, 

• otherwise z;) = 0 for all u,v G M+. 

Proof. Equation (2) implies that 

0 < Si{t) + r2S2(t) =u + r 2 V + Xi{t) + r2X2(t) + (1 - rir 2 )Li(t) + Ei{t) + r2E2{t). 

Thus if u + r 2 V + Xi{t) + r2X2{t) < 0 then necessarily Ei{t) + r2E2{t) > 0 implying that t < t. Note that 
if (9) does not hold then at least one of the Cramer-Lundberg processes Xi{t)+r2X2{t) and riXi{t) + X2{t) 
is certain to get ruined, which shows that 4 >{u, v) = 0. 

In the following assume that (9) holds, which guarantees that at least one of Hi is positive. Without loss 
of generality we assume that /ri > 0. By considering the case that the first company does not need any help 
we arrive at the bound 

^(m, 0) > P(Vf > 0 : M + Xi(f) + X2(<)r2 > 0). 

Here 2L2W i® running infimum of the process X2. It is well known in the queueing literature that 
a.s. Xi{t)/t -G Hi and X 2 (t)/f —^ min(^2,0) as f — >■ 00 . Now we see from hi + i’2H2 > 0 that a.s. Xi{t) + 
K. 2 {t )^2 —t 00 and hence the infimum of this process is finite. This yields (j>{oo,0) = 1. 

For any e > 0 choose u so large that 0) > 1 — e/2, and consider the first passage time T = infjf > 
0 : Xi{t) + X_i{t) > u} of the reflected process above level m > 0. In order to guarantee that ^(0, u) > 1 — e 
we need to choose v so large that P(X 2 (T) + 2fi{T)ri < —v) < e/2, i.e. the second company has enough 
capital for the first to reach high capital. This is clearly possible, because T < 00 a.s. □ 






5.1. From survival function to measure 

Since 4>{u, v) G [0,1] is continuous and non-decreasing in both u, ?; > 0 it defines a finite (probability) 
measure (j){du,dv) on where 4>(u,v) = Iye[o v] dy). By Fubini we have for si,S 2 > 0 the 

following relation 

F{si,S2):= [ e-"i“-*"Xdu,dz;) 

= S1S2 / ls^>u,y>v}e~''^'"~‘‘^^dxdy(f{du,dv) 

Jml JR2 

= S1S2 / e~^^^~'^^'‘^(j>{x,y)dxdy 
JRl 

= SiS2F{si,S2). 

Similarly, 

pOO 

Fi{s):= / e-^Xd7/,0) = 5Fi(5), 

Jo- 

pOO 

F2{s):= e-^X0,du) = sF2(s) 

Jo- 

for s > 0. We will refer to (j){du, 0) and d?;) as boundary measures. Note that one can easily rewrite 
the results of Proposition 1 and Corollary 1 in terms of measure transforms F, Fi,F 2 . Finally, observe that 

lim F(si,S 2 ) = ^ 2 ( 52 ), lim F(si, S 2 ) = A(si), hm F(si, S 2 ) ='/'(0,0) 

Si—>-00 S2—>-00 Si—>-00,S2—>-C>0 

and hence given F we can easily find Fi,F 2 and (/)(0,0). 


6. Independent driving processes 

Throughout this section we assume that the claim streams are independent, and thus we can focus on 
the kernel equation in Corollary 1. For this case we provide an explicit solution for the bivariate Laplace 
transform F(si,S 2 ) of the survival probability 4>(u,v). 

Even though the kernel equation is quite different from the kernel equation of a coupled processor model, 
the method of analysis from [7], building upon [11], can still be used. In the following we present a self- 
contained (apart from a few technical properties which can be found in [7]) application of this method to 
our risk problem. 

We assume that Xi and X 2 are two independent Cramer-Lundberg processes with exponents ipi[s) and 
1 ^ 2 ( 3 ), and so 'ip{si,S 2 ) = ^/’i(si) + ^’ 2 ( 52 ), see Remark 1. It is well known that '0i(s) = q for q > 0 has a 
unique positive solution, call it $^(( 7 ). Moreover, <1)^(0) = limq^o‘hi( 9 ) is a solution of i’iis) = 0, which is 
positive when fii < 0 and $i( 0 ) = 0 when fii > 0 . 

It turns out that there exists a unique non-decreasing Levy process Yi (descending ladder time process [19, 
p. 170] corresponding to Xi) such that 

i^liq) := logEe-^^-(i) = yt - (10) 

(here and in the sequel, the positive and negative parts of some entity c are denoted by c+ := max(c, 0 ) and 
c~ := — min(c, 0)). In our case Yi is a compound Poisson process (CPP) with no deterministic drift and 
jumps distributed as inf{t > 0 : Xi{t) < 0}, see also [7]. Define two two-sided CPP processes Xl and Xr 
and two constants pr and pn as follows: 

XL{t) = Y^{r^t)-Y2{t), XR{t)=Y^{t)-Y2{r2t), 

Pl = IJ-t + , PR = pt +^ 
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where we assume that Yi and Y 2 are independent. Note that pl,Pr > 0 according to the net profit 
condition (9). Observe that the corresponding Laplace exponents are given by 


M-w) ■■= logEe-“^-(i) =PL- (12) 

li’R^-w) := logEe-“'^«(i) = pr - rn T V 

where w G iR,w ^ 0. 

Finally, we introduce the Wiener-Hopf factors for Xr,pr: 

«'+(ic) = 3fJ(w) > 0, 'i’liw) = s)fj(y;) < q, 


where Cp^ denotes an independent exponential random variable of rate pr and X r,, Xp are the running 
supremum and infimum oi Xr respectively. It is well known that 




PL 

Pl - ’4)L[-wy 


w G iK, 


(14) 


see e.g. [19, Thm. 6.16]. In a similar way we define the Wiener-Hopf factors for Xr,pr. We are now ready 
to formulate our main result. 


Theorem 1. For independent Xi and X 2 satisfying (9) and Si > $^(0) it holds that 


{lfl{si) + 'lp2{s2))F{si,S2) 


P'r tp2{s2) - y2{r2Sl) 

S 2 - r 2 Si V’l(si) + 'f’ 2 {r 2 Sl) 4']^ (-!/'! (si)) 

Pl Msi)-'>PiiriS2)'^R{-if2{s2)) 
Si - riS2 1p2{s2) +i’i{riS2) 4'))(-V'2(s2)) ’ 


where 


p'l = riPi +pt - rir2P2 > p'r = ’>’ 2^2 + pf - • 

We note that the cases S 2 = r 2 Si and si = riS 2 should be understood in the limiting sense. Importantly, 
numerical evaluation of the survival probability function (j){u, v) using Theorem 1 is feasible. We refer 
an interested reader to [13, 14] for a discussion of an efficient and accurate numerical inversion of a two- 
dimensional Laplace transform and computation of the Wiener-Hopf factors, respectively. 

Proof of Theorem 1. We rewrite (8) using measure transforms from Section 5.1: 

(si - riS2)(s2 - r2Sl)(^/’l(si) -I- i>2{s2))F{si, S2) 

= (si - riS 2 ){y 2 {s 2 ) - V' 2 (t‘ 2 Si))Yi(si)s 2 
+ (S2 - r 2 Si){'f’i{si) - '0l(riS2))F2(s2)si. 

By analytic continuation and continuity on the imaginary axis we see that this equation holds for all si, S 2 
with 3?(si),3 ?(s2 ) > 0. It is shown in [7] that 4>i(s) can be analytically continued to 3?(s) > 0. Moreover, 
in this domain 3fi(<i)i(s)) > 0 for 3?(s) > 0, 4)i(s) ^ 0 for s 0, and the identity 'tlji{^i{s)) = s is preserved. 
Thus we can plug in the above equation si = 4'i(r(;), S 2 = <I> 2 (—w) for w G iR to obtain the following: 

0 = (<I>i(w) - ri$2(-w))(-r(; - y2ir2^i{w)))Fi{^i{w))^2i-w) 

+ (4’2(—w) — r2^i{w))(w — '0i(ri$2 (—w)))A(‘ 1’2(—w))4>i(ri;). 
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Assume for a moment that w 0 and multiply both sides by which using (12) leads to 


{PL - iPl{-w))- 

{PR - iPr{-w)) 


—w — '(IJ 2 {r 2 ^i{w)) 

$i(w) 

^2{-w) 


'Fi{^i{w)) = 


F2{^2i-W)). 


Using Wiener-Hopf factorization we arrive at 


PL 

PR 


w 


■ V'2(?’2^’l(w)) 


$i(w) 

{w) —w + ’tpi{ri^ 2 {—w)) 


^rH 


$2(—R') 


A($i(w)) = 

F2{^2{-w)). 


(15) 


Consider (15) and note that the Ihs is analytic in Ift(w) > 0, the rhs is analytic in ^(w) < 0, and both 
are continuous and coincide on the imaginary axis; for w = 0 this is checked by taking the limit in the 
respective half plane, see also (16) below. Hence this equation defines an entire function. Let us show 
that it is bounded and hence is a constant, call it C, by Liouville’s theorem. According to [7] the ratios 
of W-H factors are bounded in their respective half planes. The transforms Fi of boundary measures are 
bounded by 1. Finally we consider {w + 4’2{L2^iiw)))/^i{w) for ^{w) > 0; the corresponding term on 
the rhs can be analyzed in the same way. Boundedness follows from the following simple observations: 
|<i>i(i(;)| —)■ oo,tjj 2 iw)/w —)■ C 2 ,'!/l’i($i(w))/$i('u;) = w/^i{w) —>■ c\ as |ui| —>■ oo. 

Both sides of (15) are equal to some constant C, which can be identified by taking the limit w —>■ 0. 
Suppose pi > 0 then $i(0) = 0 and hence we have 


lim 

|ii;|4.0,3?(Lti)>0 


w + 'tlJ2(r2^i(w)) 
$l(w) 


t/’i(0) + r2'02(O) = Ml + ’'2M2- 


(16) 


This results in 


C = Pl{pi + i'2M2)«f>(oo, 0), 


and for /i 2 > 0 we similarly have 


C = Pr{p 2 + rifii)(j){0, oo). 


Recall that under the net profit condition (9), we have ^(oo,0) = (j){0,oo) = 1 and at least one of pi,p 2 
should be positive. Now it is not hard to check that C/pl = p'ji, C/pr = p'^. 

Finally, pick Si > <i>i(0) which implies $i('0i(si)) = Si. By considering the left part of (15) for w = 
i/>i(si) > 0 and its right part for w = —ip 2 [s 2 ) < 0 we obtain 


F^[s^) 


Pl(V'i(si) + i/’2(r2Si)) 4'+(V'i(si)) ’ 


F,[s2) = .. 

Pr(V'2(s 2) + V'l(nS2)) i^{-'tp2{s2)) 

Combining these calculations with Corollary 1 we get the result. 

Corollary 2. For independent Xi and X 2 satisfying (9) it holds that 


□ 


<^( 0 , 0 ) 


Pr 

Cl + r 2 C 2 'I'^(oo) 


p'l ^r(-oo) 

Cl + r 2 C 2 'i>2{-oo) 
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Proof. First observe that = F{X^{epj^) = 0) and similar observations hold true for the other W-H 

factors. From Theorem 1 we see that 

F(s s )= _^ 

^ S2 - r2Sl Vl/’l('Sl) + ^2(f’2Sl) ^l(Sl) + V'2(S2)/ ^';^(V’l(si)) 

P'lSI /_S 2 _ 52 A ^fl(-V^ 2 (s 2 )) 

Sl-riS2 \'f’2is2) +'fpl{riS2) V'l(5l) +'02(52)/ (“02 (S 2 )) ' 

To get the first statement, let S 2 —t oo and then si —)■ oo. The second is obtained by reversing the order of 
limit operations. □ 


6.1. Example: contingent surplus note 

Let us provide a check for a simple particular case, which admits a direct analysis. Suppose that 
^ 2 ( 1 ) = C 20 i.e. the second company is, in fact, a surplus note, which earns a constant premium and 
provides a cover for the first company. We also assume that the net profit condition (9) holds. 

Note that our system survives if and only if the one dimensional Cramer-Lundberg process Z[t) := 
Xi{t) + {p, 2 /ri)t survives when started at u + v/ri. Hence, with 4>z{x) the survival probability for the Z 
process when started at x, 


F{si,S 2 )= [ e ^^'“4'z(.u +v/ri)dudv 


poo poo 

Jo Jo 

- r - e-^^'^n^z{u)du 

52 - 5 i/ri Jq 

Pi + P'211"! / 1 _ 1 \ 

52-5i/ri V0z(5l) i>z{riS 2 ) J 
P2 + Piri f I ri 


Sl-riS 2 V0i(^ 152) + M252 ^101 (Si ) +/i 2 Sl 


(17) 


where ^pz is the Laplace exponent of Z. 

Let us check this result against Theorem 1. Note that P 2 = C 2 > 0, 02 ( 52 ) = P 2 S 2 , ^ 2 ( 9 ) = <llP 2 and so 
ip 2 (q) — 0 implying that y 2 is a zero process. Thus dif{w) = 'f'^(w) = 1, 

T+(w) = = Ee’^r(»)’'iepz, = -= - EE. —^ 

PL-rii’{[w) fj.2 + riw/^i[w) 

and similarly = prJ{ r 2 P 2 + rc/<i>i(w)). Hence for si,S 2 large enough we have 

^l(0i(5i)) _ PL r2P2Sl +0l(5i) 

^fl(0l(5l)) Pfl/^ 2 Sl +ri0i(si)’ 

According to Theorem 1 we can write 


(0l(5l) +Ai252)A(si,S2) 


P2P'r PL 

P 2 S 1 +ri'ipi{si) PR 


ripi + P2 0i(si) - 0i(riS2) 
5l-riS2 M 2 S 2 + 0l(?’lS2) 


which indeed coincides with (17), because p'rPl/pr = ?’iMi + P 2 - The latter is established by checking the 
two cases /ri > 0 and p.i < 0 separately. 
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6.2. Example: rir 2 = 1 

Assume that rir 2 = 1 which in particular holds when there are no transaction costs, i.e. ri = r 2 = 1. 
The net profit condition reduces to a single inequality + r 2 p ,2 > 0. From (2) it follows that 

S'i(t) + r 2 S 2 {t) = u + r 2 V + Xi{t) + r 2 X 2 {t) + Ei{t) + r 2 E 2 {t). 

Observe that the ruin time r is the first jump of Ei{t)+r 2 E 2 {t). Hence we have u+r 2 V+ Xi{t) +r 2 X 2 {T) < 0 
and u+r 2 V+Xi{t)+r 2 X 2 {t) > 0 for all t < t. Therefore, r is the ruin time of the classical Cramer-Lundberg 
process Z{t) = Xi{t) + r 2 X 2 {t) started in u + r 2 V = u + vjri. Similar to (17) we get 

^"(si,S 2 )= [ +v/ri)dudv 

J ]R_I_ 

_ Ml + pi-2/ri /_1 \ 

S2-Sl/?’l ^^(si) 'il’z{riS2) J 

_ Ml + ^2M2 /_1_1_ 

S2-Sir2 VV'l(si) + V’2(?'2Sl) V’l(»'lS2) + V'2(S2) 

Let us now check this result against Theorem 1. Firstly, note that pr = r 2 PL and XR(t) = Ai(r 2 t), which 
implies that the Wiener-Hopf factors for both processes are the same: = di^{w) and'I'£(i(;) = 

Hence we have 

(V'l(si) + i’ 2 {s 2 ))F{si,S 2 ) 

^ p'r ^2{s 2) - i)2{r2Si) ^ Ml -'il)i{riS2) 

S 2 - r 2 Si V'l(si) + V' 2 (?' 2 Sl) Si - ?’lS 2 '02(s2) + ■*/'l(^lS 2 ) 

^ p'r •iAi(si) + V>2(S2) Ml v^i(si) + V’2(S2) _ Mfl _ Ml 

S2 - r2Si +'!/j2(r2Si) Si - riS2 i/’2(s2) + V^i(riS2) S2-r2Si si-riS2' 

Note that p'j^ = r 2 PR = Mi + ''’ 2 M 2 and so we indeed obtain (18). 



6.3. Numerical experiments 

In this section we provide a simple numerical illustration of our main result. Theorem 1, by computing 

F{si,S2) = SiS2F(Si,S 2) = E())(esi, 6^ J, 


i.e. the survival probability for initial capitals which are two independent exponential random 

variables of rates si and S 2 respectively. The main difficulty lies in evaluating the Wiener-Hopf factors. This 
task can be accomplished using various approaches, see e.g. [14] presenting an algorithm based on the Spitzer 
identity and Laplace inversion technique. In this study, however, we use a Monte Carlo simulation algorithm 
to obtain Wiener-Hopf factors, which is relatively straightforward to implement since the processes Xr and 
Xr are CPPs. These two processes are constructed from the inverse local time processes Yi and I 2 , which 
we discuss in the following. 

Consider any of Yi, Y 2 and for the moment drop the index to simplify notation. From the general theory 
of local times, see [19, Ch. 6.2], we know that Y is a CPP with jump distribution Po(tq" € dx\TQ < 00 ), 
where tq = inf{t > 0 : X{t) < 0}. Let us determine the jump rate Ay of Y assuming that V’'(0) = m > 0: 


(q) = Ay(E(e '^'^0 [tq < 00 ) — 1) = Ay 



d>{q)p-q ^ Ay /_ 

’^$(g)(c-M) c-MV 




where we used the exit identity [19, (8.6)] and the fact that IY('^^(0) = IY(°)(0) = 1/c in that formula. 
Comparing to (10) we see that Ay = c — m- This allows to simulate the process Y by drawing from 
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Po(Tjj” S da:|r(^ < oo), which in turn reqnires to simnlate the original process and to reject the paths with 
‘long survival’. 

In our experiment we take deterministic claim sizes Ci = C 2 = 1, and put ci = C 2 = 1, Ai = 0.5, A 2 = 0.9 
and ri = r 2 = 1.1, which ensures that /ri,/r 2 > 0 and yields = s — Ai(l — e“®). First we simulate 

N = 10,000 realizations of Afi(ep^),Xi^(ep^),X^(ep^),X^(ep^) and then assign values to the Wiener-Hopf 

factors in an obvious way, e.g. see Figure 2. Now we can use Theorem 1 


Figure 2: Simulated Wiener-Hopf factors 'l'J('in),(—lu),m) (from top to bottom). 



to plot E(/)(esj, 632 ), see Figure 3 where we fix S 2 = 1 and let Si G [1,10]. Note that the survival probability 
is a decreasing function of si and thus it is increasing as a function of the mean initial capital Ecs^ = 1 /si. 
Finally, we report cj){0,0), the survival probability from zero initial capitals. Corollary 2 yields 0.281 and 

Figure 3: Probability of survival E(j!)(es;^, 6 ^ 3 ) for exponential initial capitals as a function of si with S 2 = 1. Dots represent 
direct simulation of the survival probability not using Theorem 1. 



0.277 using positive and negative Wiener-Hopf factors respectively; the direct simulation gives 0.279. 

7. Making some capital transfers impossible 

One can modify our risk model so that only the second company is allowed to help the first one, but not 
the opposite way. In other words, we regard the second company as an insurer of the first (against deficit) 
with its own stream of claims. Note that this new model is governed by the equations 

Si{t) = u Xi (t) LI(t) El (t), 

S2(t) = V + X2(t) — riLi(t) + E2{t). 
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In the following we specialize our results to this case by taking limits as r 2 —> oo. The survival probability 
we thus obtain corresponds to survival of both insurer and reinsurer. 

Firstly, the kernel equation becomes 


// N rn ^ ^ , '^{si,S2)-i’{riS2,S2) ^ , 

tp{sl,S2)F[Si,S2) = C2Fi[Si) + [ 82 ), 

Si - riS2 

which is immediate from Proposition 1 and (1). Secondly, the net profit condition now reads 

/i2 >0, /i2 + ’"iMi > 0 ; ( 19 ) 

which is seen by repeating the steps of the proof of Proposition 2. Assuming independence of claim streams, 
we obtain the following corollary of Theorem 1. 

Corollary 3. Assume that r 2 = oo and Xi,X 2 are independent and satisfy (19). Then for Si > $i(0) it 
holds that 


(t/'l(si) + 1p2is2))F{si,S2) = (ai 2 - J'lAti 


Si 


+ ^ 2 ( 1*2 + rifii) 


S2 '0l(si) - l/’l(riS2) 


Si - riS2 '02(S2) + l/’l(riS2) V'2(s2)^'i (-l/’2(s2)) ' 


Proof. Let us examine the limiting quantities (as r 2 —> 00 ) from the statement of Theorem 1. Firstly, 
stay the same, but 


'^r{w) -t 1, 




H2^2{-w) 

w 


To see this notice that pr ^ 00 and so XR{epj^) < Fi(ep^) —>■ 0 a.s., which yields the first limit. The second 
is then obtained from (14) (written for ‘R’) and (13). Finally, observe that pf = ri/ii+/^ 2 ,p)j /?'2 —>■ P- 2 —rip.f 
and take the limit in the equation of Theorem 1. □ 

Finally, let us assume that neither company can help the other, i.e. ri = r 2 = 00 . Thus we retrieve the 
standard bivariate ‘or’ problem, where the ruin means ruin in at least one line, see [8, Eq. (1.5)]. The kernel 
equation becomes 

l/>(si, S2 )E(si, S2) = C2Fi[si) + CiF2{s2). 

Now assume that fii, fi 2 > 0 and that Xi,X 2 are independent, and take the limit as ri —>■ 00 in Corollary 3. 
Noticing that 


T+(u;) ^ 


Mi^i(^^) 

w 


'SlM 1- 


we obtain 

(V’l(Sl) +'02(S2))-F'(S1,S2) 

This yields the correct product formula 


M 1 M 2 MiA^2 
V'l(si)) '02(S2)' 


E(si,S2) 


Ml T2 

V’l(si)) ■!/’2 (s2)) 


providing yet another check. 
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8. Open problems 


Various interesting directions for future work exist with respect to the present model. Firstly, one may 
try to extend the model to allow for more general driving processes, that is for Levy processes with negative 


jumps. This seems to be a hard problem contrary to many other settings, where such generalizations require 
little additional effort. In fact, there are essential difficulties in each step leading to Theorem 1: 

• Definition of the model. In general one can not use the recursive definition of Section 2.1. The definition 
in Section 2.2 (or a variant of it) seems to be a natural one. Nevertheless one has to establish that 
there is a unique solution to this modified Skorokhod problem. 

• The kernel equation. We believe that Proposition 1 holds for a general bivariate Levy process (Vi, X 2 ) 
with negative jumps. One may try to employ the generator of the underlying Markov process instead 
of an o{h) reasoning. There are various problems on this way including the required smoothness of 
the survival function. Alternatively, one may try to use a martingale approach, but the choice of an 
appropriate martingale is far from trivial. 

• Identification of boundary measures. One may follow the ideas from the proof of Theorem 1. In the 
Levy case, a highly problematic step is to bound both sides of (15) by a constant. 

Secondly, it is important to understand if the kernel equation in Proposition 1 characterizes the survival 
probability in some sense. In the case of a reflected Brownian motion in a quadrant the basic adjoint 
relation characterizes the stationary distribution together with the boundary measures, see [23, Thm. 3.4]. 
This property allowed Dai and Harrison [12] to construct an algorithm computing the stationary distribution 
based on the basic adjoint relation. 

Finally, one may try to interpret Theorem 1 (and the quite similar Theorem 1 in [7]) to provide a 
probabilistic approach. Moreover, one may consider some alternative risk measures of the model in (2) such 
as discounted external injections of capital. One may also try to extend the model to a setting with multiple 
companies. It seems that a probabilistic approach could be very helpful in this respect. 
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Appendix 


Proof of (5). Using monotonicity of (j) in both arguments we see that the rhs of (5) provides a lower bound 
for 4){u,v). The upper bound follows in a similar way. 


(ffu, v) < o(h) + (I — \h)(f){u + Ci/i, V + C 2 h) 




(f)(u + cih — X r 2 {v — y C 2 / 1 ), C 2 /i)/r(da;, dy) 



According to (4) the difference of the bounds is of order o{h), which completes the proof. 


□ 
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